Can Hyperfine Excitation explain the Observed Oscillation-Puzzle of Nuclear Orbital 

Electron Capture of Hydrogen-like Ions? 



ON 

o 
o 

(N 



o 

=5 
S3 



> 

(N 
(N 

o 
o 



Nicolas Winckler, 1 Katarzyna Siegieh-Iwaniuk, 2 Fritz Bosch, 1 Hans Geissel, 1 Yuri Litvinov, 3, and Zygmunt Patyk 2 

Gesellschaft fur Schwerionenforschung (GSI), Planckstrasse 1, D-6^291 Darmstadt, Germany 
2 Soltan Institute for Nuclear Studies, Hoza 69, PL-00-681 Warsaw, Poland 
s 'Max- Planck- Institut fur Kernphysik, Saupfercheckweg 1, D-69117 Heidelberg, Germany 

Modulated in time orbital electron capture (EC) decays have been observed recently in stored 
H-like 140 Pr 58 + and 142 Pm 60+ ions. Although, the experimental results are extensively discussed in 
literature, a firm interpretation has still to be established. Periodic transitions between the hyperfine 
states could possible lead to the observed effect. Both selected nuclides decay to stable daughter 
nuclei via allowed Gamow- Teller transitions. Due to the conservation of total angular momentum, 
the allowed EC decay can only proceed from the hyperfine ground state of parent ions. In this work 
we argue that periodic transitions to the excited hyperfine state (sterile) in respect to the allowed 
EC decay ground state cannot explain the observed decay pattern. 

PACS numbers: 23.40.-s, 31.30.Gs, 32.10.Fn 



Experiments studying electron capture decay of stored 
highly-charged 140 Pm and 142 Pm ions have been per- 
formed recently at GSI, Darmstadt. In these experi- 
ments, the highly-charged radioactive ions were produced 
via fragmentation of about 500 MeV/u 152 Sm projectiles. 
The mono-isotopic beams in a selected ionic charge-state 
were separated in-flight by the fragment separator FRS 
[l| and stored in the storage-cooler ring ESR Q. In the 
ESR, the ions were cooled employing first the stochas- 
tic 0] and then the electron cooling Q. The cooled 
ions were circulating in the ESR with revolution frequen- 
cies of about 2 MHz. Their decay properties have been 
measured with time-resolved Schottky Mass Spectrome- 
try (SMS) [H, @, 0]. Each ion at each revolution induced 
a mirror charge on a pair of capacitive pick-up plates 
installed inside the ring aperture. The subsequent fast 
Fourier transform yielded revolution frequency spectra. 
While the frequencies in such spectra reflect the mass- 
over-charge ratios of the ions 0, [3, the areas of the 
frequency peaks are proportional to the corresponding 
number of stored ions and to the square of their atomic 
charge @, H|. 

Two methods have been developed for the decay stud- 
ies [9]. In the first method, the number of stored parent 
and daughter ions has been monitored in time. A few ten 
to a few hundred parent 140 Pm and 142 Pm ions have been 
stored as fully-ionized, hydrogen- (H-like) or helium-like 
(He-like) ions 13, HI • Constants of the three-body (3 + , 
and orbital electron capture (EC), \ec, decays as 
well as the atomic-loss constant, \i oss , have been deter- 
mined, where the latter is mainly due to collisions with 
the rest gas atoms and atomic capture of the electrons in 
the cooler. In Table I we summarize the decay constants 
obtained in Refs. [3, 11 1. In the following we will refer 



12] 



to this method as the many-ion (MI) spectroscopy. 

The second method has been applied in Ref. 
and employs the sensitivity of SMS to single stored ions 
[H, @, 0] • At maximum three parent ions were simultane- 



ously injected into the ESR and their frequencies, that 
is the masses, were monitored in time. In EC decay the 
mass of the ion changes by the corresponding Q ec value 
while the atomic charge-state is preserved. Therefore, 
the EC decays were unambiguously identified by observ- 
ing sudden changes in the revolution frequency of the 
ions. Several thousands EC decays have been measured 
in this way. We refer to this method as the single-ion (SI) 
spectroscopy. Surprisingly, it has been observed that the 
number of EC decays per time unit deviates from the 
expected exponential decay. The data were described by 
adding a modulation term (1 + 0.2 x cos(0.885 • At)), 
where t is the time after the injection of the ions into the 
ESR, superimposed on the exponential decay [Hj]. The 
fit parameters taken from Ref. 12| are given in Table II. 



The interpretation of the observation above has at- 
tracted the attention of many physicists working in differ- 
ent fields. Quantum beat phenomenon due to the emitted 
neutrinos, which are flavour - but not mass eigenstates of 
the weak-interaction hamiltonian, has been suggested by 



TABLE I: Measured /3 + , EC and atomic- loss decay constants 
for fully_4onized, H-like (in bold), and He-like 140 Pr [Hj and 
142 Pm fffj ions obtained in MI measurements (see text). The 
decay constants are given in the rest frame of the ions. The 
total decay constant A = \ec + + Moss is presented in 
the last column. 



Ion 




Xec [s- 1 ] 


Moss [s ] 


A [s- 1 ] 


140 p 59+ 


0.00158(8) 




0.0003(1) 


0.0019(1) 


140 pr 58+ 


0.0016(1) 


0.00219(6) 


0.0003(1) 


0.0041(1) 


140 p 57+ 


0.0015(1) 


0.00147(7) 


0.0003(1) 


0.0033(1) 


142p m 61+ 


0.0123(7) 




0.0004(1) 


0.0127(7) 


142 pm «0 + 


0.0126(3) 


0.0051(1) 


0.0004(1) 


0.0181(3) 


142p m 59+ 


0.0139(6) 


0.0036(1) 


0.0004(1) 


0.0179(6) 
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TABLE II: The decay parameters obtained in SI measure- 
ments (see text) for 140 p r 58 + (upper part) and 142 Pm 60+ 
(lower part) [Isj]. In case a), the fits of the data points have 
been done assuming the pure exponential decay N' EC (t) — 
N(0) ■ Xec ■ e _A * (the quantity N(0) • Xec is constant), where 
A = Xec + X@+ + Xi oss and the prime index denotes the time 
derivative. In case b), a modulated in time EC decay constant 
has been assumed N' EC {t) = N(0) ■ Xec ■ e — A '(l + a x cos(cot-\- 
cp)). The decay constants X for 140 Pr have large uncertainties 
due to a short time of the total observation [12J. 



Fit parameters of 140 p r 58 + data 


Method 


A [S- 1 ] 


a 


co [s- 1 ] 


a) 


0.0014(10) 






b) 


0.0015(10) 


0.18(3) 


0.89(1) 


Fit parameters of 142 Pm 60+ data 


Method 


A [S- 1 ] 


a 


co [s- 1 ] 


a) 


0.024(4) 






b) 


0.022(4) 


0.23(4) 


0.89(3) 



several authors, see for instance Refs. [13|, [lj]. However, 
this possibility is strongly disputed, see for example Refs. 
liL 16j |. An alternative explanation has been proposed 
based on the coupling of the electron and nuclear spins 
to the rotation in the ring (Thomas precession) [T7I ] . 

In this Paper we investigate the hypothesis of periodic 
transfers between the hyperfine states (e.g. due to in- 
teractions with electromagnetic fields or with the cooling 
system) of H-like 14 °p r 58+ an( } i42p m 60+ iong ^ Due to 

conservation of the total angular momentum, the upper 
hyperfine state does not decay by allowed EC decay [10( 
which can cause the observed modulations. 

A detailed theoretical description of the EC decay rates 
in H-like and He-like ions has been performed in Refs. 
Til 19j |. There, it has been proven that the ratio of these 



rates equal to A^ likc « 3/2 • A^7 likc to a few percent. 
Based on this fact, we will show that the measured de- 
cay constants in Refs. (ToL [llT | (see Table I) cannot acco- 
modate the hypothesis of periodic transfers between the 
hyperfine states. 

Moreover, we discuss the possibility of observing the 
modulated decays also in MI experiments, which, if suc- 
cessful, would be of a great advantage due to the small 
statistics presently achievable in SI measurements. 

Both investigated nuclei 140 Pr and 142 Pm have spin 
I = l + and decay by allowed Gamow- Teller transition 
to the ground state of stable daughter nuclei with spin 
+ . The parent H-like ions can therefore have two hyper- 
fine states |— ) and |+) each with total angular momenta 
F~ = I — 1/2 and F+ = 7+1/2, respectively. The 
order of the hyperfine states depends on the sign of the 
corresponding nuclear magnetic moment. For the mag- 
netic moment parallel to the nuclear spin the hyperfine 
ground state has spin F~ and in the opposite case it is 



TABLE III: Hyperfine splitting parameters. The ion and the 
nuclear transition are given in the first two columns. Nuclear 
magnetic moments used in this work are given in the second 
column. The hyperfine splitting SE and the relaxation decay 
constants Xhf are given in the third and fourth columns, re- 
spectively. The excitation constant b (see text) is given in the 
last column. 



Ion 


Transition 


fl/flN 


SE [eV] 


A,,/ is- 1 ] 


b [s- 1 ] 


140p r 


1+^0+ 


+2.5 


1.26 


38.2 


7.8 


142 Pm 


1+^0+ 


+2.5 


1.12 


26.2 


5.4 



F + . However, in the allowed Gamow- Teller EC transi- 
tions I — » J± 1, due to conservation of the total angular 
momentum, only states can decay that have spins F . 

The energy splitting 5E between the two states |— ) 
and |+) can be estimated by using (20j : 



SE = -a(aZ) 3 — ——A(aZ)mc 2 , 



UN m p 
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(1) 



where 



/i, /ijV) an d a are the electron mass, proton 



mass, nuclear magnetic moment, nuclear magneton, and 
the fine-structure constant, respectively. The relativistic 
factor A{aZ) is defined as follows 



A(aZ) 



(2y/l -{aZf- 1)0 -(azy 



(2) 



We note, that the corrections for the nuclear charge dis- 
tribution, the nuclear magnetization distribution and the 
QED effects are not included in Eq. ([I]). 

The probability for a spontaneous transition from the 
excited hyperfine state |±) with angular momentum F^ 
to the ground state |+) is expressed by a simple for- 
mula EL El 



1 2FT + 1 



1 



(3) 



The calculated energy splitting SE and the decay proba- 
bility Xhf are presented for H-like 140 Pr and 142 Pm ions 
in Table IIIII The magnetic moment of 140 Pr has been 
estimated in Ref. j23[. The 142 Pm nucleus has a similar 
structure, therefore we applied the same magnetic mo- 
ment values as in the case of 140 Pr. Note, the hyperfine 
decay constant Xhf depends on the ninth power of the 
atomic number Z. In our case, the ground state of H-like 
ions |— ) has the spin F~ and can decay via EC and (3 + 
decay. However, the excited state |+) can only decay by 
(3 + decay. We note, that ions in both hyperfine states 
can be lost due to atomic interactions. We will denote 
the ions as active or non- active to distinguish the ions 
which can or cannot decay by the allowed EC decay, re- 
spectively. The number of active ions (A) in the state 
I — ) is NA(t) and the number of non-active (N) nuclei in 
the state |+) is N N (t) 
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We simulate the time-modulation of the number of EC 
decaying ions by introducing a simple mechanism: nuclei 
in the A-state |— ) are periodically excited with a proba- 
bility bx {l+cos(wi+<ys)} to the N-state |+). The excited 
state |+) decays spontaneously to the ground state |— ) 
with decay constant Xhf (See Fig. [T]). Experimentally 
such periodic excitations could be due to motion in the 
electromagnetic fields of the ESR or due to spin-flip re- 
actions in the cooler. The relevant hyperfine splitting 
parameters are summarized in Table IIIII The modula- 
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FIG. 1: Fractions of EC active mother nuclei (upper solid), 
non-active mother nuclei (lower solid) and the total number of 
mother nuclei (dot) and fraction of daughter nuclei (dashed- 
dotted) as a function of laboratory time. Note, the total num- 
bers of mother and daughter nuclei almost do not oscillate in 
time. 

tion of the EC decay constant is described by the set of 
four parameters {Xhf, b, w, ip}. For the sake of simplicity 
we put ip — 0. We can write three coupled differential 
equations for the number of parent A-ions N A (t), N-ions 
Npf(t), and for the number of daughter nuclei Nr>(i): 

N' A (t) = -(Xec + X loss + X 0+ )N A (t) 

+ X hf N N (t) — bx{l + cos(ujt)}N A (t), (4) 

N' N (t) = -(X loss + X p+ )N N (t) 

- X hf N N (t) + b x {I + cos(Lut)}N A (t), (5) 

N' D (t) = -X loss N D (t) + X EC N A {t) 7 (6) 

where the prime index denotes the time derivative. These 
equations have been solved numerically using the Euler 
method. A time step of 0.00056 s has been chosen. In the 
beginning we establish an approximate relation between 
the constants Xhf and b. Let us define 



< N N (t) > 



< (N N (t) + N A (t)) > 



(7) 



which is the ratio of the average number of N-nuclei to 
the total number of parent ions. An average number of 
ions decaying from the N-state to the A-state (N — > A) 
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FIG. 2: The ratio b/Xhf (dotted), the renormalization factor 
Xec{o)/Xec{0) (solid) and the average occupancy of the ex- 
cited hyperfine state (dashed) as a function of the oscillation 
amplitude a. (The meaning of the vertical axis is given in the 
legend.) 



approximately equals to xXhf and it should be equal to 
the number of ions 6(1 — x) being excited from the A-state 
to the N-state (A -» N) 



xX h f ~ b(l - x) 



blX h f 
l + b/Xhf' 



(8) 



For a fixed ratio x, the parameters A^ / and b are free. 
Using Eq. [7] the ratio x can be connected to the mod- 
ulation amplitude a: x w a/(l + a). Comparing this 
relation with Eq. El we find that a « b/Xhf For the 
observed a w 0.2 l_2j we obtain the ratio x « 0.17. 
The numerical solutions of Eqs. [3E1 ar e plotted in Fig. 
Q] We observe that, if we introduce the modulation pa- 
rameters derived from SI experiments, we are not able 
to describe the decay curves measured with MI method 
(for the fit parameters see Tables IHITT|) . The reason for 
this is simple: effectively only the 1 — x fraction of par- 
ent nuclei is active. Thus, the effective decay constant 
A w (1 - x)X EC + A/3 + A i oss is smaller since some of 
the ions are in "sterile" N-state. Therefore, if the peri- 
odic transfer between the hyperfine states occurs in the 
MI experiments, the decay constants obtained in H-like 
ions shall be readjusted. Due to the absence of hyper- 
fine splitting, no such adjustment is needed in the case 
of He-like ions. The measured ratio of decay constants 
\H-Uke i \He-hke t(J be correcte d approximately by 

1/(1 — x) = 1.20. The numerical solutions of Eqs. 0][6] 
fitted to the experimental data [l(| gives the correction 
factor of 1.23. However, the readjusted experimental ra- 
tios for 142 Pm and 140 Pr ions are 1.77(7) and 1.83(10) 
which disagree, respectively, by about 4 and 3 standard 
deviations to the theoretical ratio of 1.5. 

The numerical solutions of Eqs. [3KS1 for a fixed value 
of the oscillation amplitude a were fitted to the MI ex- 
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perimental mother and daughter populations [10]. As a 
result we obtained the ratio b/Xhf, the average occupancy 
of the excited hypcrfine state x and the EC decay con- 
stant \Ec( a )- The results as a function of the amplitude 
a are plotted in Fig. O 

Let us analyze a different model and assume that in- 
deed the electron capture probability varies in time as 
1 + 0.2 cos(0.89 At + ip), where At is the time interval 
from the creation until the decay of the ion. The origin 
of such modulations is still intensively discussed. Exper- 
iments with implanted neutral 142 Pm and 180 Re atoms, 
performed in Berkeley and in Garching [2||, respec- 
tively, have not observed time-modulations of EC decays. 

The mathematical model with the time-modulated EC 
constant Xec has the following form: 

dN M (t) = -\ E c{l + a cos(ojt + <p)}N M {t)dt 

- (Xioss + Xp+)N M (t)dt, (9) 
dN D {t) = X E c{l + acos{ut + ip)}N M {t)dt 

- X loss N D (t)dt, . (10) 

The decay constants can be taken from Table UJ and the 
corresponding time-modulated parameters a and u> can 
be taken from Table [TTJ The solution of Eq. [9] can be 
written as: 

The Eq. [10] has been integrated analytically and also 
checked numerically. The obtained results show that the 
assumption of a time-dependent decay probability ob- 
served in SI measurements is in perfect agreement with 
the decay and growth curves extracted from the experi- 
mental data of MI experiments. 

However, we can construct a different approach and 
assume that all decay constants are time-modulated or, 
what is mathematically equivalent, that the time scale 
fluctuates. The numbers of mother Njst(t) and daughter 
No(t) nuclei are connected by the two following equa- 
tions: 



We have shown that if the modulation appears only 



for H-like ions then the ratio A 



dN M (t) 
dN D (t) 



— dt{l + a cosiujt + if)} x 
(Xec + Xio SS + X fj+ )N M (t), (11) 
dt{l + a cos(ujt + ip)} x 
(\ E cN M (t) - X loss N D (t)). (12) 



Eqs ITTirrSI can be easily solved by introducing a fluctu- 
ating time scale t' connected with time t by the relation 
t' = t — a sin(ujt + ip)/uj. Thus, the numbers of mother 
and EC daughter ions can be expressed by analytic func- 
tions of time t' (see e.g. [l(| ) ■ For large t (compared 
with I /to) the ratio of both times is approaching one. 
We see that experimental populations in MI can be fit- 
ted with almost the same parameters for a ^ and for 
a = 0. 



H-lik 
EC 



;e / \ He 
l A EC 



EC Uke « 3/2 
Such an increased 



191 ] has to be multiplied by 1.23 
value is very improbable from the view of the present ex- 
perimental and theoretical data. The Af^^/Af^^ 6 
ratios from MI experiments arc in excellent agreement 
with the factor 3/2. This disagreement could be resolved 
if the described modulation appears in He-like ions as 
well. However, there are no hyperfine states in He-like 
ions and the discussed mechanism is physically disabled. 
We emphasize, that in order to describe consistently the 
present data of MI and SI experiments, the modulation 
phenomenon with similar amplitude has also to occur in 
He-like ions. Therefore, it is indispensable to perform SI 
experiments on He-like 142 Pm 59+ or 140 p r 57 + ions. Fur- 
thermore, the ratio X" ~ like / X^~ like turned out to be 
a sensitive probe and more accurate measurements are, 
therefore, required. 

We have investigated the consistency of the measured 
data in respect to the modulated EC decay constant. 
We have shown that the 20 % modulation amplitude 
determined from SI experiments translates into a tiny, 
much less than a percent, modulation amplitude of decay 
curves measured with the MI spectroscopy. Such accu- 
racy is presently out of reach for the MI experiments. 
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